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The  analysis  reported  herein  was  performed  as  an  in-house  effort  by  the 
author  for  the  Aerospace  Power  Division,  Air  Force  Aero  Propulsion  Labora¬ 
tory,  Wright- Patterson  Air  Force  Base,  Ohio.  The  work  was  accomplished 
under  Task  31402,  Project  3141,  Electric  Propulsion  Technology,  The  diagnostics 
program  was  conducted  from  June  1967  to  January  1969.  This  report  was 
submitted  by  the  author  in  February  1969. 
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ABSTRACT 


Ae  part  of  an  in-house  colloid  diagnostic  program,  the  first  detailed 
mapping  of  propellant  mass  flow  distribution  from  a  colloid  source  has  been 
completed,  and  the  program  has  shown  the  Influence  of  anisotroptc  or  non- 
uniform  flow  conditions  within  the  beam  on  the  efficiency  of  source  operation. 
The  operating  conditions  of  high  propellant  feed  pressure,  low  voltage,  and 
high  resistivity  of  the  sodium  todide/glycerol  propellant,  which  created  a  very 
low  specific  charge  beam,  caused  a  nonuniform  distribution  of  specific  charge 
and  beam  current  as  a  function  of  beam  divergence  half  angle.  The  assumption 
of  uniform  beam  conditions  would  result  in  a  12%  error  in  computing  thrust 
and  a  25%  error  in  determining  overall  efficiency. 
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SECTION  I 

INTRODUCTION 

For  about  10  years,  the  Air  Force  and  the  NatiorM  Aeronautics  and 
Space  Administration  (NASA)  have  been  conducting  research  and  development, 
both  in-house  and  contractually,  on  many  types  of  electric  propulsion  thrustor 
concepts.  The  purpose  of  this  work  is  to  develop  space  thrustors  that  will 
provide  significant  performance  gains  over  present  and  future  chemical  thrustors 
by  minimizing  sys'-em  weight  and  providing  mission  flexibility  of  thrust  vectoring, 
throttling,  precisely  controlled  impulse  bits,  and  rapid  cycling. 


Electric  engines,  unlike  the  conventional  rocket  engines,  require  an 
external  power  source  to  accelerate  propellant  gases  to  high  velocities.  This 
additional  power  requirement  and  its  associated  added  weight  are  power  and 
weight  penalties  that  must  be  offset  by  other  elect/ ic  thrustor  performance 
characteristics  or  electric  engines  will  not  be  considered  for  space  applications. 


It  is  important,  in  order  to  define  accurately  performance  gains  for 
electric  thrustors,  to  understand  thoroughly  the  characteristics  of  the  engine 
exhaust.  The  propellant  mass  flow  distribution  and  velocity  vector  must  be 
known  precisely  before  the  thrust  vector  can  be  computed  and  overall  engine 
efficiency  can  be  determined. 


The  purpose  of  this  program  was  to  determine  the  effects  of  beam  current 
and  charge-to-mass  (specific  charge)  distributions  on  propellant  mass  flow 
distribution  in  a  colloid  ion  engine  source.  These  effects,  if  any,  can  be  seen 
in  the  calculation  of  thrust,  efficiency,  and  specific  impulse. 


1 


AFAPL-TR-69-13 


SECTION  n 

THE  SPECIFIC  IMPULSE  PARAMETER 


One  of  the  main  performance  advantages  of  many  types  of  electric  engines 
is  the  ability  to  operate  over  a  broad  range  of  “specific  impulse.’’  Let  ue 
define  this  term  and  briefly  show  how  it  affects  engine  performance  (Reference 
1). 


Specific  impulse  is  the  length  of  time  in  seconds  an  engine  can  operate  at 
a  certain  thrust  level  on  a  given  amount  of  fuel.  The  higher  the  specific  impulse 
value  for  a  given  thrustor,  the  more  economical  that  thrustor  is  in  using  its 
propellant. 

Specific  impulse  is  defined  as: 

isp  =  F/Wp  ( i  ) 

where 

Isp  =  specific  impulse,  in  sec 
F  =  thrust,  in  lbs 

Wp  =  propellant  weight  flow  in  lbs/sec 
Total  system  weight  for  an  electric  thrustor  is  given  by 

w  =  p  a  +  w  t  <  2 ) 

t  t  p  p  b 

where 

total  input  power  to  the  electric  thrustor  in  kw 

power  supply  specific  weight,  in  lbs  of  power  supply  per  kw 
output 

total  system  weight,  in  lbs 
thrust  time,  in  sec 
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Since  electric  Input  power  can  be  expressed  as 

Pt  =  F  lap  ( I02/ 45.977)  (3) 

where 


77  =  electrical  efficiency  in  percent,  Equation  (2)  can  be  expanded 
and  we  obtain 

Wf  =  F I$p  *  (lOS!ap/45  917)  +  F(  tb/ lap ) 

=  Flap  *  (I02ap/45.9T7)  +  Imp/I*p  (4) 


where 


Imp  =  mission-determined  total  impulse  in  lbs-sec. 


Equations  (2),  (3),  and  (4)  show  that,  for  a  mission-determined  thrust  and 
a  specific  weight  for  a  given  power  supply,  the  total  system  weight  is  influenced 
by  the  engine  efficiency,  total  impulse,  and  specific  impulse.  Total  Impulse  is 
determined  by  the  requirements  for  a  specific  mission.  Equation  (4)  shows 
that  the  highest  engine  efficiency  is  always  desirable  in  terms  of  minimizing 
system  weight,  but  a  trade-off  occurs  between  power  Bupply  weight  and  total 
propellant  weight  as  a  function  of  Isp: 

Pf  =  flltpl 
=  f  { l/Isp) 

Figure  1  shows  a  typical  system  weight  trade-off. 


A  broad  range  of  specific  impulse  is,  therefore,  of  interest  to  the  systems 
designer  because  of  the  tradeoffs  he  can  make  in  achieving  optimum  system 
power  and  system  weight.  In  general,  these  tradeoffs  involve  the  parameters 
of  thrust  level,  specific  impulse,  power  supply  weight,  engine  efficiency,  and 
mission  duration. 
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SECTION  HI 

THE  COLLOID  ION  ENGINE 

The  colloid  ion  engine  is  one  type  of  electric  propulsion  device  that  offers 
high  potential  for  operating  efficiently  and  performing  many  useful  missions 
in  space.  The  following  paragraphs  briefly  describe  the  principle  of  operation. 

A  conducting  liquid  passing  through  a  fine  metallic  capillary  tube  is 
subjected  to  an  intense  electrostatic  field  as  it  reaches  the  tip  of  the  capillary. 
This  field  is  produced  by  raising  the  potential  of  the  capillary  with  respect  to 
a  metallic  plate  containing  a  hole,  the  center  of  which  lies  on  the  axis  of  the 
capillary.  If  the  electric  field  is  high  enough,  the  electric  forces  on  the  liquid 
exceeds  the  surface  tension  forces,  and  the  liquid  surface  ruptures  and  forms 
small  charged  droplets  called  colloids.  These  charged  droplets  are  then 
accelerated  to  high  velocities  by  the  field  and  pass  through  the  aperture  in 
the  dovmstream  electrode. 

Beyond  the  downstream  electrode,  the  charged  droplets  must  l>e  neutralized 
to  prevent  space  charge  buildup  in  the  beam  and  resultant  thrust  reduction. 
Figure  2  shows  a  simplified  schematic  of  a  colloid  engine. 

The  velocity  of  a  charged  particle  is  related  by 


w 


V/m 


=  OqIsp,  in  metar*  /fc 


where 


(5) 


V  =  applied  potential  in  volts 

q/m  =  represents  the  charge-to-mass  ratio  of  the  charged  droplet 
in  couiombs/kg 

2 

gQ  =  gravitational  constant  =9.81  meters/sec 


The  specific  Impulse  obtainable  by  a  colloid  device,  tnerefore,  depends 
not  only  upon  the  applied  potential,  but  upon  the  charge-to-mass  distribution 
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Figure  2.  Simplified  Schematic  of  Colloid  Engine 

in  the  beam.  It  will  be  shown  later  that  the  efficiency  of  a  colloid  engine,  in 
large  part,  is  also  influenced  by  this  distribution  (Equation  23). 

Figure  3  shows  the  charge- to- mass  ratio  (c)  (specific  charge)  as  a  function 
of  applied  potential  with  specific  impulse  as  a  parameter. 
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SECTION  IV 

DETERMINATION  OF  COLLOID  CHARGE -TO-MASS 
DISTRIBUTION  (ISOTROPIC) 

The  charge-to-mass  ratio  of  a  colloid  beam  is  not  uniform  from  colloid 
particle  to  colloid  particle  but  varies  from  several  coulombs  per  kg  to  many 
thousands  of  coul/kg.  A  technique,  therefore,  must  be  used  to  determine  the 
average  or  mean  charge- to-mass  ratio  of  the  beam,  c. 

1.  THE  TIME-OF- FLIGHT  SPECTROMETER 

In  1963,  H.  Shelton  and  E.  Cohen  developed  a  time-of-fllght  (TOF)  spectrom¬ 
eter  which  permits  the  specific  charge  distribution  (c)  of  the  entire  beam  to 
be  measured  instantaneously  (Reference  2).  The  device  works  in  the  following 
way  (Reference  3). 

A  shielded  collector  is  located  at  a  distance  L  from  the  colloid  source 
as  shown  in  Figure  2.  Because  of  the  electric  field  configuration  at  the  needle 
tip,  charged  droplets  reach  their  terminal  velocities  within  a  few  millimeters 
and  then  travel  at  that  velocity  to  the  downstream  collector.  This  particle 
velocity  is  given  by  Equation  (5).  Since  all  the  droplets  fall  through  the  same 
potential,  there  will  be  a  velocity  distribution  corresponding  to  the  specific 
charge  distribution  in  the  beam.  The  current  read  at  the  collector  is  the  sum 
of  each  of  the  individual  specific  charge  species. 


When  the  potential  on  the  needle  is  pulsed  to  ground,  no  more  charged 
droplets  can  be  formed,  and  a  current  decay  can  be  observed  at  the  collector. 
Particle  species  having  the  highest  value  of  c  and,  therefore,  the  highest 
velocity,  are  the  first  to  stop  contributing  to  the  collector  current.  When  the 
last  particle  of  this  specie  arrives  at  the  collector,  the  current  will  drop  to  a 
lower  level.  The  current  will  then  continue  to  decay  to  lower  levels  as  different 
c  species  arrive  at  the  collector  (see  Figure  4).  In  practice,  specific  charge 
distributions  are  equally  continuous  so  that  the  decay  is  a  continuous  function, 
rather  than  the  step  function  as  shown. 
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Figure  4.  Oecilloscope  Display  of  TOF  Decay 


The  time  at  which  the  current  decays  to  a  new  level  represents  the  time 
required  for  the  last  particle  of  that  specie  to  travel  from  needle  to  collector 
and  is  given  by  (Reference  3). 


72c V 


2  Vt 


(6) 


(7) 


This  represents  the  specific  charge  of  a  given  specie  in  the  beam. 
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By  substituting  in  Equation  (11).  we  find 


00 


171  /5v  f0  ,dI 

C  *  fCD  t 

J  t  d  1 

•  L.°  [/<>]' 


1/2 


2V 


(12) 


(13) 


Equations  (10),  (12),  and  (13)  can  be  simplified  by  observing  the  relation¬ 
ships 


and 


r2 


,  CD 


=  f  f*f ( f ) df 
0 


(  14) 


.  00 


T  =  J  t  tit  ) dt 
0 


( IS) 


where 


I  dl  r® 

^  -  -t/7dI 

X)  o 


(16) 


and 


.  oo 


1  ~  T ”  /  ,dl 
l0  0 


(17) 


We  use  these  two  relationships  (Eqs.  16  and  17),  and  Equations  (10),  (12), 
and  (13)  become 


2  V  T  2 


(10) 


1/2 


1/2 


y*7  7* 

L*  (7)* 


2  V  l  t  *  )"* 


(12) 


(  13  ) 
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If  the  beam  is  assumed  to  be  isotropic  in  relation  to  specific  charge 
distribution  and  mase  flow  distribution,  by  the  derivation  in  Reference  (6)  the 
thrust,  specific  impulse,  and  efficiency  become 

f  =  yrv  V'2  <20*1 

C  D 

(1  +  cos/9)2 

where  the  factor  17  D  (the  isotropic  divergence  efficiency  =  4  )  (24) 

arises  from  the  divergence  of  the  beam. 


_  %/Tv  1/2  1/2 

Isp  -  —  c  ’’o 


■ 


(cl'V 


The  separability  of  17  and  77^  is  a  consequence  of  the  assumption  of 
isotropy.  However,  the  worst  case  would  be  if  the  beam  were  confined  to  the 
surface  of  a  cone.  The  two  efficiency  factors  would  again  be  separable,  and 
the  divergence  efficiency  would  be 


77-.  -  cosZ/3 


(24  ) 


—  T  _  1/2 

From  the  TOF  collector,  (t)  Rnd  (t  )  may  be  determined.  Then  c,  c  , 

1/0"  ^ 

and  (c  )  may  be  calculated.  These  values  are  then  used  to  calculate  the 
performance  characteristics  of  a  colloid  engine  as  given  by  Equations  (20), 
(21),  (22),  and  (23). 


2.  DETERMINATION  OF  t  AND  t  FROM  TOF  MEASUREMENTS 

When  a  TOF  measurement  is  recorded,  the  needle  potential  is  shorted  to 
ground,  and  an  oscilloscope  records  the  current  decay  as  a  function  of  time. 
A  Polaroid  camera  attached  to  the  scope  records  the  TOF  traces.  The  TOF 
trace  may  be  divided  into  incremental  segments  along  the  time  axis  as  shown 
in  Figure  5. 


14 


A  FAPL-TR-69-13 


Equations  (18)  and  (19)  then  may  be  expanded  as  follows: 


'ft  i  So*(  &tl'  . 

"" B  L„— 


a  f  Tf  1 

t  -2  J  j. 

0  A0  nw=  1/2 


where 


Yj-  ,  ♦  Yj  , 

At  =  A*  ,  I  »  =  -i-s - L  ,  n  » 


n  i  - 1  +  nj 


and  scale  conversion  factors  are 

p  sec/cm  ,  k  =  ^tomps/cm 


.a  2(A»)2A2k 

f 

f/V'+  ni)i 

t  YJ-«+Yi)l 

'  ‘o 

L 

J  =  l 

ll  2  /' 

^  2  )\ 

2(Ai>*£*k  r 

'  (. 

Vo+  ^  + 

3  /Y.+M  5  / 

IQ  l 

2  \ 

2  / 

2  V  2  /  2  V 

I  -Nt 


a  2(A  x)Hz 


[ 


N  +N  Y.  +Y. 


+  Y,  +  2Y2  +  3Y3+  •  NYn  j 


♦  £  «'«1 

N  s  I 


In  like  manner,  where  =  total  number  of  At’s  chosen  for  calculation. 


{  h: 

0  *0  -*»«/* 


At 


n*«l/2 


(25) 


(26) 


(27) 
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t  = 


A*ki 


e  *  i 


j  - 1 


rYH*vi  • 

|  A*  ki. 

[  Y0-Y, 

l  2 

>  S  10 

l  2 

Y2  +  Y3 


v.*vi 


A*k4_  f 
J0  1 


T 


+  Y,+  Y2+  Ys  + 


'n! 


A*kA  r  Y, 


Ni 


> « *  -a.  r 

ro  1 


i  'J 


( 28) 


N  =  I 


3.  ASSUMPTIONS  MADE  FOR  CONVENTIONAL  CALCULATION  OF 

ENGINE  PARAMETERS 

The  procedures  and  equations  Just  presented  have  been  used  by  the  electric 
propulsion  Industry  to  determine  diagnostic  performance  of  colloid  engines 
composed  of  single-needle  arrays  and  multi-needle  arrays.  In  determining 
specific  charge,  thrust,  specific  Impulse,  and  efficiency,  however,  the  following 
assumptions  were  made  by  the  Industry. 

a.  Charged  particles  are  accelerated  to  their  final  velocities  in 
a  very  small  distance  compared  to  the  TOF  collector  distance. 

b.  All  charged  particles  fall  through  the  same  acceleration 
potential.  The  energy  required  to  form  the  charged  droplets  is 
small  compared  to  the  fraction  of  the  total  applied  potential. 
Therefore,  almost  all  of  the  applied  voltage  goes  into  accelerating 
the  charged  particles. 

c.  The  beam  is  isotropic  in  regard  to  specific  charge  distribution 
and  mass  flow  distribution  over  some  finite  cone  of  half  angle  (3. 

Assumption  (b)  has  been  challenged  by  Burson  (Reference  3).  He  has 
found  that  the  velocity  distribution  is  due  not  only  to  specific  charge  distribution, 
but,  to  some  degree,  to  a  voltage  distribution  caused  by  losses  during  droplet 
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formation.  He  summarizes  that  the  most  reasonable  loss  meohanism  at  the 
needle  tip  Is  Joule  heating  and  conduction  along  the  needle  but  that  this  has  not 
been  confirmed  experimentally. 


The  author  has  Investigated  the  validity  of  assumption  (c)  upon  whloh 
Equations  (20'),  (2T),  and  (23‘)  have  been  baaed.  If  the  mass  flow  la  not  Iso¬ 
tropic  over  some  finite  cone  of  half  angle  P,  then  these  equations  give  Incorrect 
values  of  thrustor  performance.  The  mass  flow  at  some  value  of  half  angle  )S 
la  given  by 


where  I^g  Is  the  current  to  a  collector  situated  at  the  half  angle /S,  and  Is 

ihe  specific  charge  determined  lrom  the  current,  l£.  If  the  mass  flow  Is 
an.'sotroplc  with  P  ,  it  can  be  due  either  to  an  anisotropic  current  distribution, 
a  variance  In  c^g  determined  from  half  angle  to  half  angle,  or  both  conditions. 
Figure  6  shows  the  possible  conditions  of  maes  flow  and  specific  charge 
distribution  with  p . 


The  fallowing  describes  the  techniques  used  and  the  equations  derived  to 
determine  the  colloid  engine  performance  calculations  based  on  anisotropic 
conditions  within  the  beam.  If  all  conditions  are  isotropic,  the  equations  still 
apply. 


17 


AFAPL-TR-69-13 


18 


Figure  6.  Possible  Current,  Specific  Charge,  and  Propella*:*  Mass  Flow  Distributions 
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SECTION  V 

DETERMINATION  OF  THRUSTOR  PERFORMANCE 
CHARACTERISTICS  BASED  ON  ANISOTROPIC  CONDITIONS 

To  determine  how  m  varies  in  the  beam,  it  is  necessary  to  determine  the 
value  ofc  and  I^a  at  Individual  points  within  the  beam,  and  then  to  perform 
a  summation  of  these  data  points  to  determine  c  (the  average  for  the  entire 
beam),  rj  ,  and  Isp, 

The  following  equations  are  derived: 

For  each  data  point,  the  thrust  component  of  a  section  of  the  beam  in 
the  x-directlon  is 


to  5  rn,,  v  o  «in  9  cos  d> 
(3a,  (3a  f3a  T 


<  29) 


The  subscript  (3a  indicates,  as  shown  in  Figure  7,  each  individual 
data  point,  e.g.,  for  each  value  of  (3  from  0“  to  70°,  a  varies  in 
incremental  degrees  from  0°  to  360?  From  Figure  7, sin  9  cos  - 
cos  (3 ;  and 

l(3a  - 


P°  ' 


2  V  c 


1/2 


-(3a 


_  1/2  !  i 

where  ,  c^  are  given  by  Equation  (101)  and  (12 * ) 

l(3  a 


f  =  ^  /2  V  c"2  cos  (3 

(3at  TTpa  [3a 

Therefore,  the  total  thrust  in  the  x  -direction  contributed  by  the  beam  is 

=  If, 


Pa 

Power  in  the  beam  is 


/3a, 


-  e-  r/3a  1/2 

2  v  X  y  cos  P 

B*  ^  0a 


IN 


V  I  I 

Pa 


Pa 


(  30) 


(31) 
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Figure  7.  Geometry  of  Source  Collector 
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and  total  propellant  mass  flow  is 


M  s  I  *  I  -s 


0a 


L-  •»«  -  L.  — 

/3a  /9a  c/9a 


(32) 


Therefore,  thrust  efficiency  due  to  specific  charge  distribution  and  beam 
divergence  is 


2  MP 


IN 


Specific  Impulse  Is  given  by 


-*„f  r  @a 

2V  2.  — — -  c_  co* 

0q  c/3q  /3a 
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2  [I  (&-)}[  *IlJ 

%K'0an[  0a  & 
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crt  COS 


Isp  = 


/3a  73a 


pa 
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I  -=? 


«0 


0a  c/3a 

and  the  mass-weighted  average  value  of  specific  charge  Is 


c  = 


|> 


I  I 


I 

M 


0a 


*  (t0a,«",I0a,a+ 


also, 


_  }  m0a  c0a 

i/a  Pa 

~T^ 


y  ^3a  *  *00*1  f  (  *0a  *2  ^ 

pa  c0a  (c/3a,i 


F  (calculated  for  non-divt.gent  beam  ) 
/2V  M 


(34) 


(35) 


(35  ) 


The  anisotropic  thrust  determination  (Equation  30)  is  simply  a  summation 

2 

of  Equation  (20’)  for  t}D  =  cos  0,  the  case  where  the  beam  Is  confined  to  the 
surface  of  a  cone.  Also,  for  small  values  of  0, 


l  +  COS0 


cos  0 
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the  isotropic  determination  of  thrust  approaches  the  hollow  cone  determination. 


In  other  words,  if  P  is  small,  Equation  (20 ‘ )  with  f) 


1/2 


'( l  +  cos  P )8 


'0  V  4 

can  be  summed  at  each  TOF  data  point  to  arrive  at  approximately  the  same 
results  as  those  obtainable  from  Equation  30.  If  the  beam  is  confined  to  a 
hollow  cone  at  some  value  of  P  ,  Equation  (20 1 )  with  t;  q  =  cos^3can  be  summed 
at  each  TOF  data  point  to  arrive  at  the  same  results  as  those  obtained  from 
Equation  (30). 
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SECTION  VI 

EXPERIMENTAL  SETUP  AND  PROCEDURES 


Since  the  experimental  plan  was  to  observe  small  segments  of  the  beam 

—  ~T7Z 

and  sum  the  results  to  determine  an  average  c  and  c  ,  for  the  entire  beam, 
a  TOF  collector  had  to  be  designed.  Three  different  design  approaches  were 
considered. 

1.  SEGMENTED  HEMISPHERICAL  COLLECTOR 

A  TOF  collector  (Figure  8a)  In  the  form  of  a  hemisphere  insulated  so  that 
incremental  currents  could  be  measured  was  considered  as  a  possible  approach 
to  mapping  the  beam.  This  approach  nad  the  advantage  of  collecting  almost 
the  entire  beam.  Also,  no  movable  parts  were  included.  However,  the  fabrication 
of  individual  “patches”  on  the  inner  surface  to  collect  the  incremental  currents 
appeared  to  be  a  formidable  task.  The  electrical  circuitry,  including  leads 
from  each  “patch”  collector,  and  the  switching  panel  connecting  the  collector 
to  the  oscilloscope  appeared  too  cumbersome  and  complex.  This  approach  was 
dropped. 

2.  BOOM  COLLECTOR 

A  TOF  collector  in  the  form  of  a  small  spherical  collector  attached  to 
the  end  of  a  boom  (Figure  8b)  was  considered.  On  the  surface,  this  seemed  a 
good  approach.  The  small  collector  3uld  be  moved  in  and  out  and  in  a  circular 
arc,  pivoting  about  the  flexible  metal  vacuum  fitting.  However,  knowing  the 
precise  location  of  the  ball  collector  was  another  matter.  A  pointer  at  the 
opposite  end  of  the  boom  could  locate  a  position  on  a  movable  plate.  With  a 
fixed  distance  x  of  the  plate  from  the  vacuum  flange,  the  bail  could  be  rotated 
in  a  circular  arc  of  radius  through  the  beam.  By  moving  the  boom  in  or  out, 
a  larger  or  smaller  circle  of  radius  a^  could  be  mapped.  Figure  9  shows  the 
position  of  the  plate  as  a  function  of  x  and  z  for  various  half  angles, /3  .  As  can 
be  seen,  many  of  the  rings  on  the  plate  are  close  to  each  other,  making  the 
pointer  very  difficult  to  position.  In  addition,  as  f3  is  changed,  x  changes,  and 
therefore,  the  TOF  distance  varies.  This  makes  the  summation  computations 
more  difficult.  Because  of  the  difficulties  this  approach  was  not  tried. 
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Figure  9.  Positions  of  Collector  Boom  Pointer  and  Movable  Plate  for 
Locating  Boom  Collector 


AFAPL-TR-69-13 


3.  GEARED  BALL  CURRENT  COLLECTOR  (BC) 

The  third  design  approach  was  to  position  a  small  (9/16"  dia.)  stainless 
steel  ball  on  the  end  of  a  gear  rack  and  pinion.  The  rack  and  pinion  were 
fastened  to  the  outside  of  a  3/4"  stainless  steel  tube  about  2l"  In  length. 
Inside  the  tube,  running  along  the  axis  of  the  tube,  was  a  1/4"  rod  which  was 
attached  to  the  pinion.  Figure  10  shows  the  assembly.  By  turning  the  rod,  the 
pinion  rotated  and  raised  or  lowered  the  ball  on  the  end  of  the  rack.  In  addition 
to  this  movement,  the  mechanism  could  be  rotated  about  the  axis  passing 
through  the  hollow  tube  by  turning  a  collector  support  brace  attachea  to  the 
tube  (Figure  11).  The  ball  collector,  therefore,  could  be  rotated  about  the 
colloid  needle  axis  a  full  360°,  aztmuthally  at  various  radii  depending  o.i  how 
far  the  rack  was  moved. 

In  addition  to  the  preceding  movements,  another  important  movement  could 
be  made.  By  turning  a  threaded  rod  which  passes  through  a  larger  (collector 
control  rod  brace)  brace  attached  to  the  vacuum  flange  and  then  passing  freely 
through  the  collector  support  brace,  the  entire  mechanism  could  be  moved 
along  the  axis,  inward  or  outward,  about  5  inches. 

With  this  arrangement,  it  was  possible  to  keep  a  fixed  distance  L  (TOF 
distance)  from  the  needles  and  position  the  ball  collector  at  any  location  in 
the  beam. 

Thus,  for  each  half  angle  oi  fi  ^  =  0,  /32-  /Jg...  70°,  the  ball  collector 

could  map  the  beam  at  azimuthal  angles  of  a  ^  =  0,  a,,,  tig,...  a  -  360°,  at  a 
constant  TOF  distance  (Figure  7).  The  ball  collector,  therefore,  could  map 
over  the  surface  of  a  sphere  L  inches  away  from  the  engine.  If  a  different  TOF 
distance  was  desired,  this  could  easily  be  changed,  and  the  collector  again 
could  be  mapped  over  a  spherical  surface. 

4.  TOTAL  CURRENT  COLLECTOR  (TC) 

In  addition  to  the  movable  ball  collector,  a  hemispherical  collector  of 
1  3/4"  radius  was  designed  and  fabricated.  When  mounted  on  a  rod  extending 
from  the  colloid  engine  vacuum  flange,  the  total  current  collector  could  be 
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Figure  10.  Geared  Ball  Current  Collector  (BC)  Assembly 
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extended  outward  from  the  flange  and  rotated  to  a  position  directly  In  front  of 
the  engine  for  a  total  current  TOF  trace.  Then  this  could  be  compared  with  the 
summation  ball  collector  (BC)  TOF  traces.  Figure  12  ahows  the  oollold  engine 
TC  collector  arrangement,  while  Figure  13  shows  the  mechanical  arrangement 
on  the  back  of  the  colloid  engine  flange  (upper  portion  of  flange)  for  extending 
and  rotating  the  TC  collector.  The  transparent  plastic  rod  shown  is  for  rotating 
the  collector  to  a  position  In  front  of  the  engine. 

5.  COLLOID  SOURCE 

A  single  needle  and  a  three-needle  array  were  used  In  the  diagnostic 
program.  The  single-needle  engine  had  a  stainless  steel  needle  of  0.008"  O.D., 
0.004"  LD.  x  4  cm.  The  needle,  during  normal  operation,  was  positioned  In  the 
aperture  of  the  extractor  electrode  along  the  electrode  axis,  and  In  the  plane  of 
the  grid.  A  mechanical  arrangement  could  move  the  extractor  electrode  slightly 
forward  and  backward  so  that  the  effects  of  the  extractor  positioning  on  the 
beam  shape  could  be  determined.  However,  for  this  study,  the  needles  were 
positioned  flush  with  the  front  surface  of  the  extractor.  Figure  13  shows  the 
mechanical  setup  for  (center  of  flange)  positioning  the  extractor.  The  extractor 
electrode  was  made  of  1/16"  stainless  steel  with  a  7/32"  diameter  aperture 
with  rounded  edges.  The  3-needle  array  was  obtained  under  contract.  The 
platinum  needles  are  0.014"  O.D.  and  0.004"  LD.  with  a  10*  funnel  level  at  the 
emitting  end. 

The  propellant  selected  for  the  tests  had  a  mixture  ratio  of  3  5  grams  of 
NAI  to  100  ml  of  glycerol.  The  purpose  of  the  sodium  iodide  is  to  produce  an 
ionlcally  conducting  liquid  which  Is  amenable  to  charging  during  the  spraying 
process  {Reference  7).  This  propellant  Is  used  widely  throughout  the  electric 
propulsion  Industry  because  of  Its  high  degree  of  solubility.  Glycerol  is  used 
because  it  has  a  low  vapor  pressure  and  Is  a  good  solvent  for  many  ion-forming 
salts.  A  low  vapor  pressure  Is  desirable  to  minimize  the  loss  of  propellant 
maBs  by  evaporation  during  particle  formation  and  acceleration  (Reference  8). 

The  propellant  feed  system  consists  of  a  fluid  reservoir  connected  to  a 
manometer  and  a  liquid  nitrogen  pressure  bottle.  Feed  flow  rate  Is  controlled 
by  a  bleed  valve  and  the  required  pressure  Is  read  on  the  manometer. 
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Figure  12.  Single-Needle  Source  and  Movable  Total  Current  Collector  (TC) 
Assembly 
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Figure  13.  Controls  for  Moving  Total  Current  Collector  Assembly  and 
<  for  Positioning  Extractor 
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The  geometrical  arrangement  of  a  single-needle  engine  and  the  TC  and 
BC  collectors  is  shown  inside  the  vacuum  chamber  in  Figure  14.  A  screen  oan 
be  seen  around  the  sides  of  the  needle.  This  screen,  which  is  eleotrloally 
connected  to  the  extractor  and  at  the  same  potential,  provides  a  barrier  pre¬ 
venting  electrons  in  the  chamber  environment  from  going  to  the  positively 
charged  needle.  The  screen  around  the  BC  collector  can  be  seen  and  is  charged 
negative  with  respect  to  the  collector  to  suppress  secondary  electrons  that 
are  emitted  as  the  oharged  droplets  strike  the  collector. 

6.  VACUUM  CHAMBER  AND  TOF  SCOPE  SYSTEMS 

The  vacuum  chamber  consisted  of  a  12"  diameter  glass  tee  and  a  6"  CVC 
vacuum  pumping  system  with  a  LN0  cold  trap.  A  reducing  flange  mated  the  12  " 

.<  ^  _g 

tee  with  the  6  diffusion  pumps.  Pressures  on  the  low  10  Torr  scale  were 
recorded  during  engine  operation. 


A  Tektronix,  Inc.  Model  355  memory  scope,  mounted  in  a  movable  power 
supply  rack,  was  used  to  record  the  TOF  traces.  The  scope,  which  was  con¬ 
nected  to  the  TOF  collector,  was  triggered  by  a  thyratron  discharge  circuit  as 
it  pulsed  the  capillary  potential  to  ground.  A  Polaroid  camera  attached  to  the 
scope  recorded  the  traces. 

Figure  15  shows  the  entire  colloid  engine  diagnostic  setup. 


On  the  left  is  the  positive  and  negative  high  voltage  power  supplies,  the  oscillo¬ 
scope/camera,  and  the  necessary  voltmeters  and  ammeters.  Through  the 
spherical  part  of  the  tee  can  be  seen  the  colloid  needle  and  extractor.  The  box 
resting  on  the  vacuum  system  controls  the  scope  resistances  and  plugs  into 
the  scope  and  either  one  of  the  TOF  collectors.  The  electrical  schematic  for 
the  entire  setup  Is  shown  in  Figure  16. 
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Figure  14.  View  of  Single-  Needle  Colloid  Source,  Total  Current  Collector, 
and  Ball  Curre  ;t  Collector,  Inside  Vacuum  Chamber 
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SECTION  VII 

RESULTS  OF  DIAGNOSTIC  TESTS 

The  first  experimental  setup  employed  a  singLe-needle  colloid  source, 
while  the  last  two  used  the  three-needle  array.  Table  1  lists  the  conditions 
under  which  the  test6  were  made. 


TABLE  I 

OPERATING  CONDITIONS  FOR  COLLOID  DIAGNOSTIC  EXPERIMENTS 


Conditions  Common 
to  Each  Test 

TEST  I 

(Single  Needle) 

TEST  n 

(Three-  Needle 
Array  ) 

TEST  HI 

(Three  Needle 
Arroy ) 

Source  Vote  rial  o nd 

Geomet  ry 
( OD  x  ID)  mils 

Stoinie  ss 
Steel 

8  x  4 

Platinum 

14  x  4  VY’th  10° 
Funnel  Level 

Platinum 
14x4  With  iQ® 
Funnel  Level 

Propellonf  (  grams  Nol/ml 
glycerol  ) 

15  /  100 

15/100 

15/100 

Propellant  Resistivity 
(  K  n  -cm  at  20 ®  ) 

2.95 

2  9 

2.9 

Needle  Voltage 
(volts) 

2500 

5300 

4800 

Extractor  voltage 
(  volts ) 

-450 

-  4„0 

-485 

Von  m  Pressure 
(  Tor  r  ) 

1.6  x  I0'6 

33  x  I0‘6 

3  x  I0'6 

Feed  System  Pressure 
(  Tor r  ) 

2.5 

2.0 

2.8 

Source  Current 
(  amperes) 

3  x  I0'6 

12  x  I0‘6 

6  x  1 0  6 

Because  of  the  great  bulk  of  data  that  had  to  lx;  reduced  and  plotted,  only 
lest  No.  II  is  reported  in  detail,  although  ample  data  is  available  for  determining 
more  performance  results  for  Test  I  and  III. 

A  digital  computer  program  was  developed  and  used  in  most  of  the  data 
reductions  and  calculations,  because  of  the  great  volume  of  data  that  needed 
to  be  reduced.  The  program  is  shown  in  the  appendix. 
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1.  TAKING  THE  DATA 

A  single  TOF  was  recorded  on  the  total  current  collector  (TC)  while  many 
traces  were  recorded  on  the  ball  current  collector  (BC).  In  taking  the  BC  data, 
a  value  of  was  first  selected  (Figure  7).  Then  the  ball  collector  was  rotated 
an  incremental  Ad  and  a  TOF  trace  was  taken. 

This  procedure  was  then  followed  through  equal  incremental  Aa ’s  over 
the  entire  360°.  The  BC  was  then  positioned  at  another  value  of  J3,  and  the 
entire  procedure  was  repeated.  TOF  traces  were  taken  at  incremental  values 
of  A  a  =  20°,  at  the  values  of  /3  greater  than  25°,  and  at  slightly  larger  values 
of  A  a  at  smaller  values  of  j3 .  This  was  done  to  prevent  overlapping  of  TOF 
readings  and  to  reduce  the  number  of  TOF  traces  that  had  to  be  reduced. 
Current  readings,  however,  were  taken  at  incremental  values  of  Aa  =  5”.  On 
all  three  tests,  approximately  500  current  readings  were  recorded;  about  270 
TOF  traces  were  made  with  135  photographed.  Of  these,  about  90  TOF  traces 
were  inspected,  and  the  data  was  reduced  for  performance  calculations. 

A  sample  BC  TOF  trace  is  shown  in  Figure  17  for  the  single-needle  Test 
No.  I  for  @  =  25°,  a  =  150°. 

The  horizontal  scale  i ,  -2  microseconds/cm  and  the  vertical  scale  is 
1  millivolt/cm  =  1.01  x  10  6  amperes/cm.  Needle  voltage  was  2500  volts  with 
3 /j.  amperes  needle  cun  ant  recorded  at  the  ammeter.  The  charge-to-masB 
ratio  determined  for  this  data  point  was  c  =  170.95  coul/kg. 


2.  TEST  II  (THREE-NEEDLE  ARRAY) 

Figure  18  shows  the  TOF  trace  on  the  TC  for  the  three-needle  array 

(Test  11),  .hilc  Figure  19  shows  a  typical  TOF  trace  on  the  BC  at  /3  25°, 

— 

a  *  30°.  The  vertical  scale  ol  the  TC  is  10  millivolts/cm  (0.92  x  10  amperes/ 

_g 

cm),  whil  that  for  the  BC  is  2  millivolts/cm  (2.02  x  10  ampere6/cm).  TOF 

distance  was  1,95"  for  the  TC  and  5"  for  the  BC.TOF  determinations  were 

c~  =  24.2  coul/kg  for  the  TC  and  c  -  66.6  coul/kg  for  this  particular  ball  collector 

(Ja 

position. 


Figure  1 
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Figure  19.  Three-Needle  Array  Ball  Current  Collector  (BC) 

TOF  Trace:  L  =  5  "  ,  (3  =  25%  a  =  30°.  cq  =66.6 
Coulombs/Kg,  c  =  63.8  Coulombs/Kg 


a.  Determination  of  Current  and  Mass  Flow  Distribution 
and  Total  Mass  Flow  Rate 

While  many  current  and  TOF  data  points  were  taken  as  a  function  of  /3  and 
a  ,  it  ib  obvious  that  a  large  fraction  of  the  l>eam  was  not  sampled  by  the  ball 
collector,  as  shown  in  Figure  20.  It  was  necessary,  therefore,  to  plot  the 
available  data  points  and  perform  a  graphical  integration  to  obtain  the  correct 
values  of  current  and  mass  flow.  This  procedure  was  as  follows. 

Figure  21  shows  the  beam  current  measured  at  the  ball  collector  for  a 
given  value  of  /3  =  25°  as  a  function  of  a  (Aa  =  5°). 

Figure  22  shows  the  TOF  charge -to-mass  calculation  for  /3  =  25",  as  a 
function  of  a.  These  points  were  determined  for  incremental  Aa  =  20°.  Straight 
lines  were  connected  to  each  data  point.  The  resultant  propellant  mass  flow 
measured  at  the  ball  collector  is  show-n  in  Figure  23  for  /3~  25°,  as  a  function 
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Figure  21.  Beam  Current  Measured  at  Ball  Current  Collector  as  a 
Function  of  A7.inv.  thal  Angle,  a  ,  for  /3  =  25“ 


Figure  22.  Specific  Cliarge,  Determined  from  Ball  Collector 
as  a  Function  of  Azimuthal  Angle  a  for  (3  =  25° 


Figure  23.  Propellant  Mass  Flow,  Calculated  at  the  Ball  Collector  as  a 
Function  of  Azimuthal  Angle,  a  for  B  =  25° 
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of  a  (  A  a  =  5°).  The  average  mass  flow  density  measured  by  the  ball  collector 
as  a  function  of  the  angle  is  then  determined  a6  follows: 


mPt, 


_!_  y 

2ir  % 


Aa  ,  in  kg  /sec 


where 


a  =  cross-sectional  area  of  ball  collector  =  0.25  in  . 


m„  =  mass  flow  measured  at  ball  collector  in  kg/sec 
Pa  each  point. 


l0a 

— —  at 

c/3a 


The  average  mass  flow  density  at  a  given  position  in  the  beam  was  deter¬ 
mined  bv  weighting  the  m  n  vs  a  graph.  Using  the  information  that  the  graph 

(3a 

o  lit* 

paper  weight  C.06  grams/in. ,  and  the  size  of  the  graph  was  9  x  6  ,  the  weight 
of  the  graph  paper  is  determined  to  be  54  x  0,06  grams  =  3.25  grams.  There¬ 
fore,  3.25  grams  =  12  x  10-10  kg/sec  (the  value  of  the  6"  ordinate  scale)  x  2v 
radians  (the  value  of  the  9"  abscissa  scale)  .  .  1  =  23.25  x  10  10  kg- radians/ 
gram-sec.  Therefore, 


m 


p(3 


/  1  \ 

/  1  _\ 

/23.25kQ-rodions\ 

V  2ir  radioes/ 

'0  25  in2/ 

\  grom-sec  / 

where  rh^ga  Aa  is  cut  out  and  weighed  on  a  Right-A-Weight  balance  scale. 

This  procedure  was  used  for  (3  =  O',  7.5',  15°,  25°,  35°,  50°,  and  70°  and 
then  a  smooth  curve  was  drawn  in  Figure  24  for  mp^g  versus  /3.  This  curve 
now  represents  the  mass  flow  density  distribution  as  a  function  of  (3 . 

The  mass  flow  distribution  as  a  function  of  (3  and  the  total  mass  flow'  can 
now  be  calculated  as  follows: 


M  - 


A  A 

•  n 


73 


(  37) 
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in  lO*10  Kilograms  /  inch2- second 
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(3  ,  Oegrees 


Figure  24. 


Propellant  Mass  Flow  Density.  Calculated  at  the  Ball  Collector, 
as  a  Function  of  Azimuthal  Angle,  a  ,  for  /5  =  20’ 
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where  AA  =  approximate  elemental  surface  area  of  circular  strip  of  radius  It, 
and  width  AS,  at  angle  /3  (See  Figure  25). 


M 


:  J  2  irr  z  tin  )3  A  )3,  mKQ/sec 

p  03 


For  the  ball  collector,  the  angle  subtended  by  the  ball  was  A/3  =  2.5°  = 
.044  radians,  and  the  distance  of  the  needle  to  the  ball  collector  was  r  -  5"  . 


M  =  6.84  Y.  sin^  •  ,n 

(3  09 

Total  current  measured  for  Test  II  by  the  above  procedure  was  determined 

*  6 

to  be  1.4  x  10  amperes.  Total  mass  flow  determined  graphically  at  each  point 

-8 

and  summed  with  the  use  of  Equation  (32)  was  2.2  x  10  kg/sec.  However,  it 
will  be  shown  later  in  Section  VII,  Paragraph  2b,  that  the  current  and  mass  flow 
must  be  corrected  because  interference  of  the  ion  beam  with  the  extractor 
electrode  and  the  collector  screen  reduces  the  percent  of  ions  reaching  the 
collector  and  causes  secondary  electrons  to  travel  from  the  screen  to  the 
collector,  further  reducing  the  measured  current  reading.  This  however,  will 
not  change  the  mass  flow  distribution.  The  corrected  mass  flow  and  the  mass 
flow  distribution  are  shown  as  a  function  of  half  angle  /3  in  Figure  26  (See 
Section  VII,  Paragraph  2b(2),  for  corrected  mass  flow  determination). 


The  average  charge-to-mass  ratio  for  ihe  beam  is  (from  Equation  35) 


63  8  coul  f  Kg 


The  dotted  line  in  F  igure  26  is  the  mass  flow  distribution  if  one  were 
simply  to  take  a  current  distribution  and  divide  by  63.8  cou’/kg  To  determine 
the  mass  flow  distribution  in  this  way  would  Lx;  to  assume  that  c  was  constant 
at  all  positions  in  the  beam  and  could  be  determined  by  a  single  TOF  reading. 
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_  1/2 

Figure  27  shows  the  variation  of  thu  parameters  and  c p  ,  as  a  function 
of  (3  .  It  has  previously  been  shown  that  the  specific  charge  varied  widely  with 
azimuthal  angle  (see  Figure  22).  There  is  a  gradual  increase  in  the  parameters, 
with  (3  out  to  an  angle  of  18°,  then  a  gradual  decrease  out  to  /3  =  70°,  The 
results  of  Figure  27  were  used  when  the  solid  line  in  Figure  26  wa6  drawn. 
Because  no  beam  current  was  measured  at  angles  of  /3  =  0°  to  7.5°  and  only  a 
neglible  current  was  measured  at  /3  =  50°  to  70°,  c p  ,  c ■  and  the  current 
had  to  be  interpolated  in  these  regions  by  smooth  curves.  Therefore,  this  part 
of  the  m  curve  in  Figure  26  is  not  as  accurate  as  the  curve  for  (3  =  7.5°  to  [3  = 
50°.  However,  only  a  small  fraction  of  the  mass  flow  exists  outside  this  region. 
The  error  in  mass  flow  is  estimated  to  be  about  2%. 

The  correct  mass  flow  distribution  shows  a  shift  of  about  5°  from  the 
incorrect  (c  =  constant)  determination. 

As  can  be  seen  from  Figure  26,  the  beam  is  hollow  along  the  engine  axis 
with  the  flow  peaking  at  the  half  angle  (3  =  36°  and  decreasing  rapidly  to  zero 
at  70°.  Figure  28  is  a  photograph  of  a  model  of  the  mass  flow  distribution  of 
the  beam  determined  by  grapns  similar  to  the  one  shown  in  Figure  23.  The  fact 
that  the  beam  is  hollow  was  verified  by  observing  the  pattern  of  propellant 
depositing  on  the  TC  collector,  as  shown  in  Figure  29.  The  pattern  on  the 
collector  is  a  result  of  an  accumulation  of  three  individual  colloid  runs,  all 
wmi  low  rp:  uF,  cha:„  ^.ns:  however,  or.  s’l  operations,  the  beam 

was  essentially  hollow.  Although  the  photographers’  lights  caused  a  glare  at 
the  center  of  the  hemispherical  collector,  little  or  no  deposit  was  observed. 

While  the  mechanics  of  colloid  particle  extraction  are  not  completely 
understood,  we  know  that  the  detachment  point  for  particles  at  the  needle  rim 
and  the  departure  angle  depends  upon  the  field  configuration  at  the  needle 
(Reference  8,  10)  The  field  influences  the  shape  of  the  fluid  meniscus  in  the 
interior  of  the  needle  and  the  shape  and  departure  angle  of  the  fluid  streamlets 
leaving  the  meniscus.  Experimental  evidence  indicates  that  spraying  occurs  in 
a  finite  number  of  discrete  jets  from  the  meniscus,  since,  by  varying  the 
needle  voltage,  a  spray  consisting  of  several  jets  can  be  obtained.  Since  it  has 
been  observed  tliat  a  decrease  in  needle  voltage  causes  beam  spread,  the  shape 
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Figure  28.  Model  of  Three-Needle  Array  Propellant  Mass  Flow  Distribution. 
Rods  lenjths  are  Proportional  to  Propellant  Flow  Rate 
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Figure  29.  Pattern  of  Propellant  Deposited  on  Total  Current  Collector 
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of  the  beam  was  probably  influenced  by  a  varying  field  strength  and  meniscuB 
shape,  (Reference  10.) 

b.  Diagnosis  of  Beam  Characteristics 

(1)  Correlation  of  Current  Readings;  Determination  of 
Droplet  Size  and  Mas6 

The  procedure  outlined  in  Paragraph  2a  increased  the  current  reading 
determined  at  the  ball  collector  from  0,55  microamperes  to  1.4  microamperes. 
However,  since  the  needle  current  mieroammeter  read  12  microamperes,  and 
the  total  current  collector  recorder  was  4.35  microamperes,  a  large  fraction 
of  the  current  was  still  unaccounted  for. 

A  determination  was  made  of  the  screening  effect  of  the  secondary  electrode 
suppressor  screen  that  surrounded  the  BC  collector.  Two  methods  were  used 
to  determine  the  screen  transparency.  The  hemispherical  screen  was  made  of 
No.  20  mesh,  0.009"  dia.  stainless  steel.  With  this  fact  known  and  assuming 
that  at  the  distance  of  5"  from  the  source  the  streams  of  charged  particles 
passed  parallel  to  each  other  and  perpendicular  to  the  screen,  a  graphical 
determination  could  be  made  of  what  fraction  of  the  beam  passed  through  the 
grid. 


By  cutting  out  the  projected  areas  cf  the  impinging  beam  on  the  ernes 
section  of  the  collector  and  weighing  the  samples,  the  average  screen  trans¬ 
parency  was  determined  to  be  about  57%.  This  result  was  checked  later  by 
using  a  Photon  Meter  (Figure  30)  to  ascertain  the  drop  in  the  intensity  of  light 
passing  through  the  screen.  The  ratio  of  the  intensity  of  light  passing  through 
a  hollow  tube  to  the  intensity  of  light  passing  through  the  tube  with  the  wire 
screen  mounted  on  one  end  was  to  Therefore,  the  average  transparency 
of  the  screen  is  the  inverse,  55-62%,  which  agrees  well  with  the  graphic  deter¬ 
mination. 

The  total  current  collector  also  had  a  screen  in  front  of  it.  While  the  ball 
curient  collector  screen  was  convex  spherical  toward  the  impinging  beam, 
the  total  current  collector  screen  was  concave  spherical  and  geometrically 
presented  a  higher  transparency  to  the  beam. 
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However,  before  equations  could  be  set  up  relating  all  system  currents, 
it  was  necessary  to  consider  further  the  geometry  of  the  ball  collector  and 
screen  (Figure  31).  The  ball  screen  had  a  larger  cross-sectional  diameter 
than  did  the  hall  itself  (about  0.8  Inches  compared  to  0.57  inches)  and  was 
physically  connected  to  a  metal  plate  behind  the  ball  collector,  and  electrically 
Isolated  from  the  ball.  This  geometrical  and  electrical  arrangement  meant 
that  when  the  Ion  beam  reached  the  ball  collector  system,  part  of  the  Ions 
Impinged  on  the  portion  of  the  screen  outside  the  ball  diameter,  and  another 
fraction  of  the  beam  passed  through  this  part  of  the  screen  and  struck  the  back 
metal  plate. 

These  ion  currents,  outside  the  ball  collector  cross-sectional  area,  while 
not  adding  to  the  ball  collector  current,  produce  secondary  electrons  which 
go  to  the  collector.  The  presence  of  electron  current  to  the  collector  has  the 
effect  of  cancelling  out  an  equal  ion  current  to  the  collector.  Therefore,  the 
ball  collector  would  receive  a  greater  percentage  of  electrons  because  of  the 
current  striking  the  entire  screen  and  the  back  metal  plate.  It  was  also  deter¬ 
mined  that  in  Region  1  the  screen  transparency  was  about  0.69  while  in  Region 
it  was  0.42  (overall  screen  transparency  was  0.56). 

Figure  31  illustrates  what  is  happening.  The  ion  current  in  Region  II 
striking  th:  screen  at  (13)  and  (C>  produces  additional  secondary  electrons 
that  go  to  the  collector. 

Equations  relating  all  system  currents  could  now  be  set  up.  An  electrometer, 
originally  placed  in  the  extractor  circuit,  was  not  functioning  properly  and 
had  to  be  removed.  Unfortunately,  no  microammeter  had  been  placed  in  the 
BC  and  TC  circuits,  which  obviously  seriously  hampers  correlation  of  currents. 

An  estimate  of  the  currents,  however,  can  be  made  by  calculating  a  value 
for  the  secondary  electron  coefficient  y.  In  Reference  11,  y  has  been  deter¬ 
mined  experimentally  to  be  equal  to  1.25  at  5300  volts  needle  potential  and  at 
a  particle  specific  charge  of  7000  coulombs/kg.  As  »  first  approximation,  we 
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Figure  31.  Schematic  of  Bali  Current  Collector  Showing  Regions  of  Charged 
Particle  Impact 
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con  state  that  the  secondary  electron  coefficient  Is  a  function  of  the  kinetic 
energy  of  the  Impinging  ion  for  constant  voltage,  or  y  =  f  (me).  Therefore, 

| 

I 

i  where 

[  X2  =  1-25 

Cj  =  63.8  coul/kg 
c^  =  7000  coul/kg 

where  m,  the  mass  of  each  droplet,  can  be  determined  by 

«  kg  4  s  3  3 

m  s  yoVol  =  1.26  x  10  — j  x  irr  =  5.  28  x  10  r 
4'  m  ■* 

Droplet  size  can  be  estimated  with  the  techniques  described  in  Reference  5, 
where  droplet  size  is  determined  by  considering  the  suiface  energy  forces,, 
ooulombic  forces,  and  field  emission  effects.  Using  the  Raleigh  criterion  for  r 
(Reference  5)  for  lew  specific  charge  particles,  we  find  rx  -  10* '  meters  = 

looo  A 

Therefore 

m  =  5.28  x  lC''8kg 

i  Using  the  field  emission  limit  equation  for  ions  for  high  specific  charge  parti- 

f  eles,  we  find 


/me,  \ 

i  a‘ 


(39) 


Therefore, 


r2  -  3  x  10  8  meters  -  300  A 


,-ie, 


m2  -  O  l‘  ‘10  kg 
By  substituting  in  liquation  (39),  we  find 


v  =  0  42 
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Using  the  calculated  values  of  screen  transparencies,  the  secondary 


electron  coefficient  and  the  cross-sectional  area  relationship, 
following  current  relationships  can  be  found. 
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which  is  based  on  the  assumption  that  all  currents  In  Region_lLproduce  secondary 
electrons. 

By  manipulating  the  above  relationships  wo  can  find 
lor* 
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The  calculated  currents  are  shown  tn  Figure  32  and  are  listed  In  Table  II. 
It  should  be  noted  that  while  the  total  electron  current  from  the  entire  screen 
is  considered  because  all  of  it  impinges  on  the  ball  collector  only  the  ion 
current  in  the  cross- sectional  area  of  RegtonJJarea  a)  is  Included.  All  previous 
summations  to  determine  the  total  ball  collector  current  over  the  sperical 
surface  o.  1  to  (3  =  70°  were  made  as  a  function  of  the  ball  collector  cross- 
sectional  area  a,  and  not  the  screen  cross-sectional  area  b. 


By  using  the  equations  for  current  relationships,  current  continuity  can 
be  shown  to  hold  tn  the  range  of  0  £  y  Si,  that  is,  for  y  >  1  several  calculated 
current  values  become  negative.  Also,  for  any  value  of  y  <  1,  the  calculated 
values  of  the  total  ion  current  emitted  from  the  needle  differ  by  less  than  10%. 
The  final  calculated  value  of  current  determined  by  using  the  calculated  value 
of  y  -  0.42  is  I*  .  *  10.82  microamperes  (Figure  32). 


(2;  f  inal  Determination  of  Mass  Flow  Rate 
The  mass  flow  rate  determined  from  Equation  32  is  now 

•  10  82  *  lO*6  -s 

M  s  - rT?S —  =  170  *  '0 


(3)  Charge-to-Mass  Ratio 

Charge-to-mass  ratio  (specific  charge)  as  determined  on  the  ball  collector 
was  o’  =  63.8  coul./kg  at  V  =  5300  volts  and  propellant  feed  pressure  of  two 
inches  of  mercury.  Propellant  resistivity  for  the  15  grams  Nal/100  ml  glycerol 
solution  was  2.9  KX2-cm  at  SOX.  It  has  been  shown  (Reference  5  and  8)  that 
the  specific  charge  is  influenced  by  needle  voltage,  needie  geometry,  feed 
sy&tern  pressure,  and  the  propellant  resistivity.  Generally,  specific  charge 
increases  with  an  increase  in  needle  voltage,  and  a  decrease  in  feed  pressure 
and  resistivity,  which  ie  caused  by  an  increased  concentration  of  Nal  in  glycerol. 
The  influence  of  capillary  potential  on  specific  charge  for  a  constant  feed 
pressure  is  shown  in  Figure  33. 
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Figure  33.  Influence  of  Needle  Voltage  on  Specific  Charge  for  Constant 
Feed  System  Pressure  (TC,  TOF,  Trace) 
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It  was  impractical  to  operate  the  needles  at  voltages  higher  than  about  6  kv 
because  the  TOF  current  traces  on  the  ball  collector  became  negligible  or 
negative.  It  is  believed  that,  at  higher  voltages,  the  increased  secondary 
electron  current  from  the  ball  screen  to  the  ball  masked  out  the  ion  current. 

The  specific  charge  recorded  on  the  TC  collector,  at  the  same  operating 

conditions  (same  pressure)  as  for  the  mass  of  data  recorded  at  the  ball  collector, 

was  c  =  24.2  coul/kg.  The  discrepancy  between  this  value  and  the  ball  collector 

reading  (c  -  63.8  coul/kg)  is  unexplained.  However,  instantaneous  flow  rates, 

* 

as  determined  by  TOF  data,  can  vary  from  the  measured  average  by  as  much 
as  40%  (Reference  12).  A  large  number  of  time-of-flight  measurements  made 
under  constant  pressure  and  temperature  conditions  will  agree  closely  with 
the  measured  average.  Therefore,  it  iB  believed  that  the  determination  of  c  by 
use  of  the  summation  ball  collector  TOF  traces  is  more  accurate  than  a  single 
total  collector  TOF  trace.  However,  it  would  be  an  interesting  experiment  to 
determine  if  there  is  any  change  in  specific  charge  with  TOF  distance.  TC  TOF 
distance  was  1,95"  and  BC  TOF  distance  was  5". 

(4)  Neutral  Particle  Efflux 

Although  no  determination  of  neutral  particles  efflux  was  made,  long  time 
tests  conducted  by  the  electric  propulsion  industry  (Reference  11)  have  shown 
that  neutral  efflux  is  only  the  order  of  1-2%  and,  therefore,  is  considered 
negligible.  These  tests  were  performed  by  comparing  average  instantaneous 
TOF  mass  flow  rate  calculations  with  the  mass  flow  rate  determined  by  actually 
weighing  the  propellant  reservoir  before  and  after  a  long  term  operation  and 
dividing  by  the  time.  This  is  defined  as  the  “measured  flow  rate  average.’’ 

(5)  Beam  Divergence  and  Anisotropic  Distribution  Effects 
on  Needle  Performance 


Thrust  and  current  and  mass  flow  were  calculated  similarly,  i.e.,  by  the 
use  of  Equation  (30)  and  the  graphical  techniques.  True  thrust,  taking  into 
account  the  anisotropic  conditions  of  specific  charge  and  current  and  also  the 

♦measured  average  is  discussed  in  the  subparagraph  (4)  “Neutral  Particle 
Efflux.’’ 
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beam  divergence,  i*>  F  =  =  23.3  micropounds.  The  true  thrust  in  the  axial 

direction,  contributed  by  the  beam  as  a  function  of  divergence  angle,  is  shown 
in  Figure  34.  If  beam  divergence  i.  neglected,  Equation  (20)  (with~  =  63. 8, c1^2  = 
7.37)  may  be  used,  and  we  get  F  =  28.8  micropounds. 

If  the  mass  flow  for  a  thrust  of  28.8  micropounds  were  considered  to  be 
distributed  uniformly  out  to  a  half  angle  of  (3  =  38°  (see  Figure  26),  the  thrust 
would  be  26.1  micropounds. 

If  one  were  to  consider  that  specific  charge  was  isotropic  In  regard  to 
divergence  angle  /3  ,  then  thrust  for  isotropic  "c,  c1^2,  and  divergent  beam  is 
F  =  23.9  micropounds. 

•l  To 

This  points  out  that  the  anisotropic  conditions  of  c  and  c  '  as  shown  in 

Figure  27  were  insignificant  at  the  positions  in  the  beam  where  beam  current 

was  significant,  in  that  they  caused  little  change  in  calculated  thrust  (-2.5%). 

_  f/2 

That  is,  although  there  was  a  variation  in  c  and  c  as  a  function  of  beam 
angle  /3,  it  occurred  in  regions  where  the  beam  current  was  small. 

The  fact  that  the  beam  could  be  considered  as  truly  hollow  as  verified  by 
using  Equation  (20')  (where  “  =  63.8,  c1/2  -  7.3,  and  (3  =36°)  to  compute  the 
true  thrust. 

The  result  was  F  =  23.2  mlcropounde  compared  to  the  graphical  deter¬ 
mination  of  F  =  23.3  mtcropounds.  The  experimentor,  however,  would  be 
unlikely  to  guess  by  a  simple  current  distribution,  or  by  observing  a  beam 
pattern  on  the  collector  such  a  precise  angle  as  (3  =  36°.  However,  this  angle 
can  be  seen  as  the  angle  where  the  mass  flow  Is  greatest  (Figure  26). 

-  l/2 

In  conclusion,  the  condition  of  anisotropic  c  and  c  had  a  negligible 

effect  in  computing  thrust,  and  since  the  beam  was  hollow.  Equation  (30),  the 

—  l/2 

summation  thrust  equation  for  anisotropic  c  and  c  '  ,  current,  and  divergent 
beam,  can  be  replaced  by  Equation  (20<).  However,  if  the  observer  assumes 
that  the  team  is  uniformly  distributed  out  to  some  f3 ,  rather  than  being  hollow, 
his  thrust  calculation  is  considerably  off  (26.1  micropounds  compared  to  the 
true  value  of  23.3  micropounds). 
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Figure  34.  Thrust  Distribution  in  Axial  Direction  Contributed  by  Elemental 

Surface  Areas,  AA,  Situated  at  Different  Beam  Divergence  Angles, 0 
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The  true  value  of  thrust  being  known,  we  can  now  find  the  performance  of 
the  three-needle  array. 


We  can  define  overall  efficiency  for  a  thrustor  operating  In  a  space  en¬ 
vironment  where  tnere  are  no  collector  screens  to  retard  Ion  flow  and  cause 
secondary  electrons  as 

Vf  * 

V, 


IONS 


w+wv-v  * 

ext 


where 


Cons  =  10*  82  +1°"6 

Ifrom  =  1,19  x  10"6  amP6 
ext 

V  =  needle  voltage  =  5300  volt 
\  =  extractor  voltage  =  -450  volts 

7}o  --  T)  F  ( 0.9) 

where  77  p  is  given  by  Equation  (33) 


•’  %  *  09 
total  input  power  is 


lnt-1  ‘ 


49  2  % 


(40) 


P  =  I,...-  V  +  I.  V-V  )  =  64  2  *  10  3  wott. 

t  IONS  from 

•  x  t 

and  specific  impulse  (Equation  34)  is 


lap  -  — —  =  62  seconds 

*  -8 

As  found  previously,  M  =  17.0  x  10  kg/sec 


These  are  the  performance  parameters  arrived  at  by  assuming  an  aniso¬ 
tropic  divergent  beam. 
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The  comparison  of  true  engine  performance  with  performance  calculated 
by  not  considering  anisotropic  conditions  and  beam  divergence  is  shown  in 
Table  III,  Table  IV  summarizes  the  calculated  three-needle  array  performance. 
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TABLE  m 

THREE- NEEDLE  ARRAY  PERFORMANCE,  AS  DETERMINED  FROM 
BALL  CURRENT  COLLECTOR  TOF  TRACES 


Performonce  Porometers 

Results  Obtained 

Needle  Current,  In 

12  fx  amps 

Positive  Current,  I*t0,Q, 

1  0  62  ^tomps 

Meon  specific  charge.c 

63  6  coul  /  kg 

~1TS 

Mean  root  specific  charge,  c 

7  3  coul  *  kg 

Total  Input  power, 

64.2  *  10  5  watts 

Propellant  mass  flow,  M 

17  0  x  10  #  kg  /sbc 

True  thrust,  Fx 

23  3  fi  lbs 

Specific  impulse,  Isp 

62  seconds 

Thrust  efficiency, 

54  7  % 

sroll  Efficiency, 

49.2% 
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SECTION  VIII 
CONCLUSIONS 


From  the  analytical  and  experimental  work  performed,  the  following 
conclusions  can  be  made: 

1.  Under  the  stated  operating  conditions  of  fluid  resistivity,  propellant  feed 
pressure,  and  needle  voltage,  the  generation  of  a  low  specific  charge  colloid 
beam  results  in  an  anisotropic  propellant  mass  flow  distribution  that  is  in  the 
form  of  a  hollow  cone  with  a  divergence  angle  of  /0  =  36°.  This  distribution  is 
caused  primarily  by  a  conical  current  distribution  and  is  influenced  to  a  lesser 
extent  by  a  specific  charge  distribution  which  showed  ~c  increasing  from  the 
axis  of  the  beam  out  to  20°  and  then  gradually  decreasing  out  to  70°.  These 
results  show  that  assumption  (c)  in  Section  IV,  Paragraph  3,  made  by  researchers 
in  colloid  propulsion,  was  not  correct  for  these  operating  conditions,  as  specific 
charge  was  shown  to  be  anisotropic,  although  it  had  little  effect  on  engine 
performance.  Mass  flow  distribution,  which  Is  also  anisotropic,  cannot  be 
Ignored. 

2,  Under  the  above  stated  operating  conditions,  the  calculation  of  engine 
performance  can  be  made  with  Equation  (30)  by  detc  'mining  the  current  dis- 

—  1/2  _  ~  7/2 

trlbutlon  and  by  holding  c  and  c  constant.  Both  c  and  c  can  be  determined 
from  a  TOF  total  collector  positioned  at  the  same  distance  from  the  needles 
as  the  summation  ball  collector.  For  hollow  beams,  Equation  (30)  may  be  replaced 
by  Equation  (20‘)  where  rj  =  cos2  /3  . 

3,  The  anisotropic  nature  of  the  specific  charge  had  negligible  effect  upon  the 
calculation  of  engine  performance  (5%  error  in  overall  efficiency),  because 
the  greatest  deviation  of  specific  charge  from  the  average  occurred  where 
beam  current  was  a  small  fraction  of  the  total  current. 

4.  The  propellant  mass  flow  distribution,  which  was  in  the  nature  of  a  hollow 
cone,  had  a  significant  effect  upon  the  calculation  of  engine  performance. 
Assuming  a  uniform  mass  flow  distribution,  rather  than  the  correct  distribution, 
resulted  In  a  12%  error  in  thrus  1  calculation  and  25%  error  in  determining 
overall  efficiency. 
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SECTION  IX 

RECOMMENDATIONS 


The  following  recommendations  are  made: 

•  The  influence  of  propellant  feed  pressure  and  needle  voltage 
upon  the  formation  of  a  hollow  beam  should  be  investigated 
because,  from  a  performance  standpoint,  a  hollow  beam  Is  the 
least  efficient  of  beam  shapes, 

•  A  determination  should  be  made  of  the  effects,  if  any,  of  time- 
to-flight  distance  on  resultant  TOF  specific  charge  determination. 

•  While  specific  charge  increases  with  needle  voltage,  it  Is  not 
known  whether  the  voltage  increase  has  any  effect  on  the  specific 
charge  distribution  with  beam  half  angle;  although  it  has  been 
observed  that  the  current  distribution  becomes  more  collimated 
along  the  beam  axis. 

•  It  is  possible  that  at  higher  voltages  the  specific  charge  may 
still  exhibit  its  anisotropic  nature,  but  this  phenomenon  may 
now  occur  in  regions  of  high  current  where  it  could  more  signifi¬ 
cantly  influence  mass  flow  distribution. 
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APPENDIX 

DIGITAL  COMPUTER  PROGRAM  DEVELOPED  FOR 
CALCULATING  COLLOID  ENGINE  PERFORMANCE 
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»1  *UPRLEMENTARY  NOTH 

u  sponsoring  kiutary  activity 

Air  Force  Aero  Propulsion  Laboratory 

Wright- Patterson  Air  Force  Base,  Ohio 

1J.  AB»T*ACT 

As  part  of  an  tn-house  colloid  diagnostic  program,  the  first  detailed 
mapping  of  propellant  masB  flow  distribution  from  a  colloid  source  has  been 
completed,  and  the  program  has  shown  the  influence  of  anisotropic  or  non- 
uniform  flow  conditions  within  the  beam  on  the  efficiency  of  source  operation. 
The  operating  conditions  of  high  propellant  feed  pressure,  low  voltage,  qpri 
high  resistivity  of  the  sodium  todide/glycerol  propellant,  which  created  a  very 
low  specific  charge  beam,  caused  a  nonuniform  distribution  of  specific  charge 
and  buam  current  as  a  funetton  of  beam  divergence  half  angle.  The  assumption 
of  uniform  beam  conditions  would  result  in  a  12%  error  in  computing  thrust 
and  a  25%  error  in  determining  overall  efficiency. 
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